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A"^ • Abstract 

^vj I The current status of experiments on the d = 2 and d = 3 random-exchange and random- 

field Ising models, as realized in dilute anisotropic antiferromagnets, is discussed. Two areas of 
current investigation are emphasized. For d = 3, the large random field limit is being investigated 
and equilibrium critical behavior is being characterized at high magnetic concentrations. 
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1 Introduction 



r^ ' The Ising model is one of the most studied and basic models for phase transitions. In this article, 

the current status of experimental studies characterizing two classic models of second-order phase 
transitions in short-range interaction systems in the presence of quenched disorder, the random- 
exchange Ising model (REIM) and the random- field Ising model (RFIM), is presented. The discussion 
concentrates on experiments in dilute, insulating, anisotropic antiferromagnets, the systems that 
have yielded the best understood data for these two models. The REIM is realized in zero magnetic 
field and the RFIM with a field applied along the spin-ordering direction. The REIM is rather 
'■^ , well characterized experimentally, theoretically and through computer simulations. The d = 2 

C ' RFIM is fairly well characterized, although the scaling behavior of scattering near the destroyed 

phase transition is still being investigated. The understanding of the RFIM for d = 3 is not as 
complete, though significant progress has been made in the past few years, and it is this model that 
will be the main focus of this short review. The early history of the d — 3 RFIM was fraught with 
controversial interpretations of the data, a result of severe nonequilibrium effects. Nevertheless, some 
experimental groups realized from the start that underlying the observed, complicated behavior is a 
^^ ' new kind of phase transition. Efforts to characterize the new critical behavior were thwarted by the 

^^ , severe nonequilibrium effects. These nonequilibrium effects have recently been overcome by going 

On ' to sufficiently high magnetic concentration and a complete characterization of the universal d — 3 

RFIM critical behavior is possible and underway. The most recent static critical behavior will be 
compared to results from computer simulations and theory. In addition to these low-field behaviors, 
much has been learned about the high-field limit of the RFIM. An overview will be given of the 
C^ ' phase diagram and the different behaviors observed. 

Experiments have been performed on REIM and RFIM systems for more than two decades. Since 
the experiments performed some time ago have been reviewed previously ||l|, they will be included 
here only as needed to give a perspective on the current physical understanding of the models. 
Theory and simulation results will be included as needed for the interpretation of the experiments. 
^ , Another classic model of ordering in the presence of disorder, the spin glass, will be covered in 

another reviewB in this Ising Colloquium and so will not be discussed here in detail, although some 
spin-glass-like behaviors do occur at low magnetic concentrations and at high magnetic fields. 
^\^ • Table 1 shows the most frequently measured static critical behaviors associated with a phase 

M , transition. We will make reference to the universal parameters defined in Table 1 as needed. 

2 Experiments on Pure d = 2 and d = 3 Anisotropic Antifer- 
romagnets 

Observations of asymptotic static critical behavior in the pure d = 2 and d = 3 Ising antiferromagnets 
are very well documented. The magnetic specific heat (Cm) critical behavior has been characterized 
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Table 1: The asymptotic forms for commonly measured static critical behaviors. The superscript + 
(-) on the amplitudes signifies T > Tc (T < Tc). The exponent values and the amplitude ratios are 
universal quantities that depend only on the general properties of the system. 

using the optical linear birefringence techniques[|[ ^ on the d = 2 Rb2CoF4 systemg and the d = 3 
FeF2 system H. The measured pure Ising critical exponents a and the amplitude ratios A+/A~ are 
in superb agreement with very many theoretical and simulation results. The birefringence technique 
is particularly useful and more accurate than pulsed specific heat techniques since it is insensitive 
to the large phonon contributions which are particularly difficult to handle for d = 2. 

The critical behavior of the staggered susceptibility and correlation length have been determined 
with neutron scattering for K2COF40 for d = 2 and in FeF2|2] for d = 3. In general, the scattering 
line shapes for the pure and REIM systems away from the Bragg scattering point follow the scaling 
behavior of the spin-spin correlation function 
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where 7 = v(2 — if). For both d — 2 and d — 3, the scaling functions used in data analysis are 
approximate onesP, nfl] that differ significantly from the mean-field (MF) Lorentzian 
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as can be seen in Fig. 1 where various scaling functions are compared. The deviation from the 
Lorentzian is more pronounced for d = 2 and for T < Tc{H) in both dimensions. 

The order parameter critical behavior has been determined using neutron scattering [[7| in the 
d ~ 2 compounds. The Mossbauer technique |l^ was used in the study of the d = 3 system. 

The results for the pure Ising model in d = 2 and d = 3 are summarized in Table 2. Note that 
the Rushbrooke scaling relation 

2/3 + 7 + a > 2 (3) 

is satisfied as an equality for both cases. Included in Table 2 are the results from a few theoretical 
and simulation studies. No attempt is made to review the vast literature on the pure Ising models. 



3 Random-Exchange Experiments in Dilute d = 2 and d = 3 
Anisotropic Antiferromagnets 

The REIM is realized in dilute, anisotropic insulating antiferromagnets when the site dilution does 
not result in strongly frustrated bonds (that would lead to spin-glass behavior). Random- exchange 
phase transitions are observed in d = 2 and d = 3 systems and these appear to be in good accord 
with theory and simulations. 

The d — 2 REIM C„i critical behavior was observed using the birefringence technique |y] on the 
magnetically dilute antiferromagnet Rb2Coo.85MgQ ]^5F4. The approximately logarithmic divergence 
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Table 2: The pure d = 2 and d = Z Ising static critical exponents obtained from experiments, 
theory and Monte Carlo simulations. 
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Table 3: The d = 2 and rf 
simulations and theory. 



3 REIM Ising static critical exponents obtained from experiments, 



is compatible with theoretical predictions [|l7[ [l8|] . The scattering critical behavi or ||l9| of the com- 
pound Rb2CoxMg]^_xF4 was analyzed using approximate scattering line shapes|9|, |lO| developed for 
the pure d = 2 Ising model. The successful analysis using these line shapes suggests that the correct 
line shape is close to the pure one. The static critical behavior of the d = 2 Ising model is quite well 
characterized by experiments and theory as shown in Table 3. 

The d = 3 REIM is similarly well characterized with birefringence, neutron scattering and 
Mossbauer experiments employing FexZni_xF2 with the results shown in Table 3 along with some 
theoretical and simulation results with which they agree very well. The critical behavior of the 
specific heat of the d = 3 REIM, measured with birefringence techniques pi|, is shown in Fig. 2. 
Monte Carlo simulations P7[ based on the FexZni_xF2 system are shown in Fig. 3. The birefringence 
technique yields a negative specific heat exponent a as predicted ||28| , consistent with a universal- 
ity class different from the pure Ising model where a is positive. Note that, just as in the pure 
case, the birefringence technique is consistent with pulsed specific heat techniques, though the latter 
technique suffers from greater concentration gradient sensitivity p9| and the large phonon specific 
heat component. The critical behaviors of the staggered susceptibility and correlation length were 
determined from neutron scattering experiments pa. The order parameter critical behavior was 
determined from Mossbauer studies [p5[. The scattering line shape scaling functions are not known 
from theory and were therefore determined directly from the scattering data in neutron scattering 
experiments pO| using Feo.93Zno.07F2. The results shown in Fig. 1 clearly indicate that the scaling 
functions are fairly close to those of the pure d = 3 case. 

The REIM universal static critical parameters for d = 2 and d = 3 are shown in Table 3 along 
with theoretical and simulation results. In both dimensions the agreement is excellent. Note that 
the Rushbrooke scaling relation (Eq. 3) is satisfied as an equality for the REIM. 
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Table 4: The d 
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3 RFIM Ising static critical exponents obtained from experiments, simulations 



4 The d = 3 Magnetic Percolation Threshold Concentration 

As the magnetic percolation threshold concentration, Xp, is approached from above in zero field, the 
equilibrium phase transition is expected to approach zero temperature. For a particular magnetic 
structure, Xp depends on what interactions exist between the different neighboring spins. For exam- 
ple, for FexZni_xF2, Xp = 0.245 provided only the dominant interaction between the body-center 
and corner spins is considered [p4[. Except very close to Xp, the much smaller interactions can be 
ignored. Close to Xp, however, the smaller interactions may drastically affect the behavior and even 
prevent ordering above Xp when they frustrate the dominate interaction. Both the extremely slow 
dynamics near percolation p5| and the sensitivity to tiny frustrating interactions pq| can cause the 
system to exhibit spin-glass-like behavior. 
A good deal of effort has focused on the 
system has a small frustrating interaction 
in experiments by Montenegro, et al.pT 



roperties of the system FexZni_xF2 for x near Xp. This 
The spin- 
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dass-like properties were first elucidated 
M. For H > and x = Xp, there exists 
a boundary that resembles a de Almeida-Thouless boundary with curvature T — To ^ H^/'t> where 
(p = 3.4, a typical spin-glass value. It was shown with neutron scattering |44| that there is no 
antiferromagnetic long-range ordering below this boundary. Much of the behavior is very reminiscent 
of a canonical spin glass. The detailed behavior of this sample has been studied experimentally |3q, Q 
and extensively modeled in local mean- field [^6[ ^ and Monte Carlo simulations pq|. 



5 d = 2 Random- Field Behavior 



Scaling arguments for Zeeman and domain wall energies by Imry and Ma|49] as well as considerations 
by Binder |5Q] leave little doubt that the d — 2 Ising transition is destroyed by the introduction of 
arbitrarily small random fields. Birefringence [p| and neutron Bragg scattering [Bl] experiments bear 
this out; no sharp phase transition is observed in equilibrium, though the rounded transition exhibits 
the expected scaling behavior. The equilibrium region is separated from a lower temperature region 
of strong hysteresis observed in the difference between data obtained upon heating after cooling in 
zero field to low temperatures and then applying the field (ZFC) and upon simply cooling the sample 
in the field (FC). The boundary separating these regions is time-scale dependent ]5^. The domain 
dynamics induced with the application of a magnetic field as well as those remaining after the field 



is removed at low temperatures have been studied experimentally p3|, p4| and theoretically pq| . 

6 d = 3 RFIM behavior for Xp < x < Xe set low fields. 

The behavior for concentrations between Xp and the percolation threshold concentration for va- 
cancies, Xe = 1 — Xp, with a relatively small applied field occupied the bulk of early experimental 
efforts PI . Much of the controversy over interpretations of experimental data involved this region of 
concentration and fields. As a result of the equivalence |5g, ^ of the dilute anisotropic antiferro- 
magnet in small fields and the random-field ferromagnet often studied theoretically, it was believed 
that concentrations near x — 0.5 would yield strong random-field effects in reasonably small fields 
and would be the best realizations of the d — 3 RFIM for phase transition studies. At the time of 
the first experiments Ipq, p9[ , |60| , |6l[ , it was generally believed, based on many theoretical arguments, 
that no phase transition would be observed. Indeed, early neutron scattering FC experiments, by 
Yoshizawa, et al.||6^, seemed to bear this out. In particular, a resolution- limited Gaussian Bragg 
peak does not occur upon FC, although subsequent experiments |54| show that the samples retain 
long-range order below the phase boundary if ZFC. In contrast, the first Cm studies^, |6^ yielded 
compelling evidence for a fundamentally new phase transition governing the behavior. The phase 
boundary T — Tc{H) ^ H^/'t' behaves as predictedlQ with (j) = 1.42 ± 0.03 for random-exchange to 
random-field crossover [p3. 

Despite the sharp Cm peak, many have argued against it as evidence of the existence of a phase 
transition. The most recent of these is the "trompe I'oeil transition" phenomenological model [p5[. 
Among the assumptions of this model are that the birefringence and Cm experiments do not yield 
the same behavior, the uniform magnetization is reflected by the square of the staggered magne- 
tization, and conventional scaling is inoperative. This phenomenological model was shown to be 
inconsistent pql when all available data are considered. In contrast, more proven and conventional 
techniques of analyzing the experimental data, as described in this review, have been very fruitful 
in providing consistent results in a meaningful scaling context. 

It is clear that the phase transition underlying the behavior is obscured by nonequilibrium be- 
havior below the boundary, Teq{H), lying just above the phase transition and scaling in the same 
manner, albeit with a slightly larger amplitude ||67[] . The equilibrium behavior above Tf,q{H) can 
be used to extrapolate the scattering data to the obscured phase transition boundary. When done 
carefully, the boundary determined in this way coincides with that determined via Cm experiments, 
which are much less sensitive to the non-equilibrium behavior that distorts the neutron scatter- 
ing data. For a long time, the nonequilibrium experiments represented the best random-field results 
available. One of the particularly interesting predictions [Q of the change in the critical behavior in- 
duced by the random fields is that the order parameter critical exponent f3 should decrease from 0.35 
to a value near zero. This can only be measured below Tc{H), i.e. in the nonequilibrium region, so it 



was not clear what would be observed. Experiments on thin films were made[68 for x = 0.52, well 
below Xe ~ 0.76. Not surprisingly, the results were peculiar. The curvature of the Bragg intensity 
versus T was such that it would require (3 >> 0.5, which is hard to justify theoretically. Magnetic 
x-ray scattering data showed similar behavior near surfaces in bulk samples, though they were inter- 
preted under the "trompe I'oeil" phenomenology and the Bragg scattering was not separated from 
the fluctuation scattering pq]. For some years it appeared that the problems of metastability below 
Teq{H) could not be avoided, i.e. that they were intrinsic to the random-field behavior as realized 
in dilute antiferromagnets. However, insight into the origins of the metastable domains finally led 
to experiments pOJ at high magnetic concentration as a way to avoid the nonequilibriunr behavior, 
as discussed later. 

The metastable domains formed upon FC are themselves quite interesting and their dynamics 
were studied in some detail both experimentally ]69|, 70, ^ |2[ ^ and through simulations JT^, |7^ , 



7 RFIM behavior for Xp < x < x^, at high fields. 

The general phase diagram features, shown in Fig. 4 for the RFIM at high fields were investigated in 
pioneering pulsed-field magnetization measurements p6| in FexZni_xF2- Low temperature single spin 
flips and the phase boundary are shown in Fig. 5 and, interestingly, the behavior of the upper phase 
boundary appears to be different for x < Xe and x > Xe- This is consistent with the differentiation 
of the behaviors observed in neutron scattering experiments above and below Xe ■ 

In recent years, it has become clear that weak RFIM (small applied field) and strong RFIM 
regimes exist for Xp < x < Xe- The Feo.31Zno.69F2 system exhibits Q typical low-field behavior for 
H < 1.5 T, with TN-TciH) and TN-Teq{H) scaling as H'^/'l' with <j) w 1.4. At larger fields, however, 
the curvature for T^ —Teq{H) changes to ~ 3.4, a value close to that observed in spin glasses. Some 
of these features were suggested qualitatively ||7q, [79[ . While the lower region has been shown to have 
antiferromagnetic long-range order upon ZFC ||80|| , no long-range antiferromagnetic order is observed 
at the higher fields below Teq{H). Instead, spin-glass-like behavior is observed. This is clearly the 
same type of behavior observed for all fields at the percolation threshold concentration|8J. The 
same type of distinctive low and high field behaviors have been observed for x as large as 0.60. The 
two regions are separated by an equilibrium boundary [p2[ p3| , as observed for x = 0.56 (Fig. 6) 
and 0.60, and it appears to decrease towards the H = boundary at finite temperature well below 
Tc{H) and approach the phase transition line at finite field, a point separating the sharp transition 
observed at low field and the more glassy transition at higher fields. This also is consiste nt p^ with 
specific heat peaks that are very sharp at low fields and quite rounded at high fields]^, |65|. 

The distinction between high- and low-field behavior behavior is observed as well in ac suscep- 
tibility experiments. At low fields, there exists a single peak which seems to be associated with 
extremely slow dynamics, either from activated dynamics or at least power-law behavior with a very 
large dynamic exponent pq, |87| . There is little hysteresis between the ZFC and FC procedures at 
low H. At larger fields the peak splits in the ZFC procedure only, with a sharp peak at slightly 
lower temperatures than the broader peak|^, Q. This splitting appears to be associated with the 
upper region of the phase diagram corresponding to the spin- glass- like behavior [p3. These effects 
have not been investigated for x > Xp. The high field region for these concentrations is still an open 
area for research. 

8 d = 3 RFIM Equihbrium Critical Behavior for x > Xe 

It was, of course, realized very early that the metastable domain walls at low magnetic concentra- 
tions took advantage of vacancies. What was not fully appreciated was that the Imry-Ma domain 
wall energy argument p3 is not applicable when domain walls can to a great extent pass through 
vacancies, avoiding the energy cost of breaking magnetic bonds. With sufficient vacancies, i.e. for 
X < Xe, domain walls can take advantage of vacancies to such an extent that the domain wall en- 
ergy can be insignificant. Interestingly, every experiment that could have detected low temperature 



hysteresis, particularly neutron scattering and capacitance |67| experiments, was done for x < 0.72, 
which is below Xe = 0.76. Higher concentrations were avoided since the generated random fields 
are small, resulting in relatively very narrow asymptotic random-field critical regions around Tc{H). 
Nevertheless, concentrations well above Xe are necessary to study the equilibrium critical behavior 
and require high magnetic fields and very fine temperature resolution. 

Figure 2 shows the specific heat data for Feo.93Zno.07F2, measured with optical linear birefrin- 
gence. The specific heat was also measured to demonstrate that, in agreement with theory[p| and 
contrary to the so-called 'trompe I'oeil' phenomenology [pSf , the data from the both techniques yield 
the same critical parameters. An important advantage of the birefringence technique is its relative 
insensitivity to the concentration gradients that tend to smear the transition |22l. Interestingly, the 
critical behavior appears to be very similar to that of lower concentrations where metastable domains 
dominate the scattering behavior. The specific heat was studied using Monte Carlo simulations p^ 



based on the Feo.93Zno.07F2 system and the results are shown in Fig. 3. Although the simulations 
are not of sufficient resolution to extract the critical exponent, the similarities in the shapes of both 
the REIM and RFIM indicate the same qualitative change from the asymmetric cusp at _ff = to 
the nearly symmetric peaks at iJ > 0. This result, however, is not seen in all MC simulations pi[ . 

The neutron scattering experiments on Feo.93Zno.07F2 show no hysteresis below the phase tran- 
sition, in stark contrast with samples with x < Xe- There is no evidence that domains form upon 
ZFC or FC in this concentration range and the line shapes are independent of the thermal cycling 
procedure, implying equilibrium conditions. 

Neutron scattering experiments on this sample were difficult to analyze since the RFIM line 
shape is not known from theory. In general, two different scaling functions are involved with the 
form 

x,(g) = A±.j"-V(g/«) + S^«''-'g(gA) . (4) 

with two independent sets of critical behavior exponents. However, simulations |22|| and high tem- 
perature series expansions ||89| strongly suggest a simpler scenario. The predictions are that the 
exponents are simply related, fj being twice 77, and the new scaling function g{q/ n) is, to a good 
approximation, the square of /{q/n). 

Not only were the universal RFIM critical parameters obtained in this experimental study, but 
the scaling analysis yielded the spin-spin correlation scaling function f{q/n). This scaling function 
is compared to several other known spin-spin correlation scaling functions in Fig. 1. Note that the 
d = 3 RFIM one seems the furthest away from the MF Lorentzian of all the examples for T > Tc{H), 
whereas for T > Tc{H) the pure d = 2 case is further than the d = 3 RFIM, though both are very 
far from MF behavior. 

The critical parameters for the RFIM with x > x^ are shown in Table 4. Certainly more effort 
is needed to complete the experimental entries and to find reconciliation between the simulation 
and experimental results. Note that some sets of exponents from the simulations violate|33] the 
Rushbrooke scaling relation in Eq. 3. 

9 The Vacancy Percolation Threshold Concentration 

The regions of low temperature nonequilibrium behavior and equilibrium behavior have been shown 
to be separated at relatively small H by a nearly vertical sharp boundary aX x k, Xe = 0.755 in 



FexZni_xF2 using Monte Carlo studies|90 . Figure 7 shows simulations on lattices of size 2 x L^ 
with L = 64 modeled closely after FexZni_xF2. Hysteresis is observed upon FC and ZFC for 
X < Xe but not above. The hysteresis for x < Xe increases for larger lattices or slower thermal 
cycling, showing that it is not simply an artifact of the simulations not being run long enough. 
The concentration dividing equilibrium and nonequilibrium behavior is very close to or equal to 
the vacancy percolation threshold concentration x^ = 0.755. Apparently, the percolation vacancy 
structure facilitates the formation of domain walls. 

To further investigate this boundary, recent experiments have been done^l| on a sample with 
X = 0.76, just above Xe- No evidence of domains has been observed for small H. Since earlier 
experiments ||67| for x — 0.72 in FexZni_xF2 gave clear evidence for domain formation, including 
a reversal of the Bragg intensity curvature just below Tc(H) upon ZFC, the boundary must be 
0.72 < Xe < 0.76, in agreement with the MC simulations. There have been no theoretical studies 
reported explaining the existence or nature of this boundary. 

10 The Current Situation and Outlook 

The recent measurement of the equilibrium critical behavior in the random-field Ising model has 
side-stepped the great difficulties encountered in the interpretation of data below the transition that 



are obfuscated by nonequilibrium phenomena in many studies at lower concentrations. Certainly, 
there is wide agreement that a phase transition exists. Although experiments for x > Xe are much 
more difficult since the random-field region is very narrow, they are being done. Interestingly, the 
experimental results are not in agreement with much of the theory and simulation results, unlike 
the REIM and pure Ising model. A reliable characterization of some aspects of the d = 3 RFIM 
universality class behavior remains to be completed. The phase diagram of FexZni_xF2 has proven 
to be quite rich in detail. An important area of the phase diagram for which a good understanding 
is being developed is the large random field limit for Xp < x < Xe- With the the progress being 
made along these two lines of inquiry, a rather complete experimental characterization of the RFIM 
in dilute antiferromagnets seems near at hand. 

Recent work has been supported by Department of Energy Grant No. DE-FG03-87ER45324. 

References 

[1] D. P. Belanger and A. P. Young, J. Mag. Mag. Mater. 100, 272 (1991); D. P. Belanger, in 
"Spin Glasses and Random Fields" , edited by A. P. Young, (World Scientific, Singapore, 1998), 
p. 251; T. Nattermann, in "Spin Glasses and Random Fields", edited by A. P. Young, (World 
Scientific, Singapore, 1998), p. 277. 



[2 
[3 

[4; 
[t; 

[8 
[9 

[lo; 

[11 

[12 
[13 

[14 
[15 
[16 

[17; 
[18; 

[19 



P. Nordblad, Ising Centennial Colloquium, 2000. 

J. Ferre and G. A. Gehring, Rep. Prog. Phys. 47, 513 (1984). 

D. P. Belanger, A. R. King and V. Jaccarino, Phys. Rev. B 29, 2636 (1984). 

I. B. Ferreira, A. R. King, V. Jaccarino, J. L. Cardy and H. J. Guggenheim, Phys. Rev. B 28, 
5192 (1983). 

P. Nordblad, D. P. Belanger, A. R. King, V. Jaccarino, and H. Ikeda, Phys. Rev. B 28, 278 
(1983). 

R. A. Cowley, M. Hagen and D. P. Belanger, J. Phys. C 17, 3763 (1984). 

D. P. Belanger and H. Yoshizawa, Phys. Rev. B 35, 4823 (1987). 

H. B. Tarko and M. E. Fisher, Phys. Rev. B 11, 1217 (1975). 

M. E. Fisher and R. J. Burford, Phys. Rev. 156, 583 (1967). 

G. K. Wertheim and D. N. E. Buchanan, Phys. Rev. 161 478, (1967). 

L. Onsager, Phys. Rev. 65, 117 (1944). 

D. P. Belanger, P. Nordblad, A. R. King, V. Jaccarino, L. Lundgren and O. Beckman, J. 
Magn. Magn. Mater. 31-34, 1095 (1983). 

M. Campostrini, A. Pefissetto, P. Rossi and E. Vicari, |cond-mat/9905078| (1999). 

R. Guida and J. Zinn- Justin, J. Phys. A 31, 8103 (1998). 

E. Brezin, J. C. LeGuillou and J. Zinn- Justin, Phys. Lett. 47A, 285 (1974). 



E. S. Sorensen, |cond-mat/0006233| (2000). 

R. Folk, Yu. Holovatch and T. Yavors'kii, Phys. Rev. B 61, 15114 (2000). 

M. Hagen, R. A. Cowley, R. M. Nicklow and H. I. Ikeda, Phys. Rev. B 36, 401 (1987). 



[20] V. S. Dotsenko and V. S. Dotsenko, J. Phys. C 15, 495 (1983); J. Phys. C 15, L557 (1983). 

[21] Z. Slanic and D. P. Belanger, J. Magn. Magn. Mater. 186, 65 (1998). 

[22] H. G. Ballesteros, L. A. Fernandez, V. Martin-Mayor, A. M. Sudupe, G. Parisi, and J. J. Ruiz- 



Lorenzo, Phys. Rev. B 58, 2740 (1998); also see A. Pelissetto and E. Vicari, cond-mat/0002402 



(2000) and D. V. Pakhnin and A. I. Sokolov, |cond-niat/991207l| (1999). 

[23] J. Wang and D. P. Belanger, unpublished. 

[24] S. A. Newlove, J. Phys. C 16, L423 (1983). 

[25] N. Rosov, A. Kleinhammcs, P. Lidbjork, C. Hohenemser and M. Eibschutz, Phys. Rev. B 37, 
3265 (1988). 

[26] D. P. Belanger, A. R. King, and V. Jaccarino, Phys. Rev. B 34, 452 (1986). 

[27] W. C. Barber and D. P. Belanger, ICM2GGG. 

[28] A. B. Harris, J. Phys. C 7, 1671 (1974). 

[29] D.P. Belanger, A.R. King, LB. Ferreira, and V. Jaccarino, Phys. Rev. B 37, 226 (1987). 

[30] Z. Slanic , D. P. Belanger and J. A. Fernandez-Baca, Phys. Rev. Lett. 82, 426 (1999). 

[31] H. Rieger, Phys. Rev. B 52, 6659 (1995); H. Rieger a nd A. P. Young, J. Phys. A 26, 5279 
(1993); J. Machta, M. E. J. Newman and L. B. Chayes, |cond-mat/0006267| (2000). 



[32] U. Nowak, K. D. Usadel and J. Esser, Physica A 250, 1 (1998). 

[33] The only measurement is using the dilation technique at lower concentration where the system 
is not in equilibrium by C. A. Ramos, A. R. King, V. Jaccarino and S. M. Rezende, J. de Phys. 
49, C8-1241 (1988). 

[34] M. F. Sykes and J. W. Essam, Phys. Rev. 133, A310 (1964). 

[35] C. L. Henley, Phys. Rev. Lett. 54, 2030 (1985). 

[36] B. W. Southern, A. P. Young and P. Pfeuty, J. Phys. C 12, 683 (1979). 

[37] W. C. Barber and D. P. Belanger, Phys. Rev. B 61, 8960 (2000). 

[38] S. M. Rezende, F. C. Montenegro, M. D. Coutinho-Filho, C. C. Becerra and A. Paduan-Filho, 
J. de Phys. C8, 1267 (1988). 

[39] F. C. Montenegro, M. D. Coutinho-Filho and S. M. Rezende, Europhys. Lett. 8, 382 (1989). 

[40] F. C. Montenegro, J. C. O. de Jesus, F. L. A. Machado, E. Montarroyos, and S. M. Rezende, 
J. Magn. Magn. Mater. 104-107, 277 (1992). 

[41] F. C. Montenegro, U. A. Leitao, M. D. Coutinho-Filho and S. M. Rezende, J. Appl. Phys. 67, 
5243 (1990). 

[42] J. H. de Araujo, J. B. M. da Cunha, A. Vasquez, L. Amaral, J. T. Moro, F. C. Montenegro, 
S. M. Rezende and M. D. Coutinho-Filho, Rev. Brasil. de Fisica 21, 115 (1991). 

[43] F. C. Montenegro, S. M. Rezende and M. D. Coutinho-Filho, J. Appl. Phys. 63, 3755 (1988). 

[44] D. P. Belanger and H. Yoshizawa Phys. Rev. B 47, 5051 (1993). 



10 



[45] K. Jonason , C. Djurberg , P. Nordblad and D. P. Belanger Phys. Rev. B 56, 5404 (1997). 

[46] E. P. Raposo, M. D. Coutinho-Filho and F. C. Montenegro, Europhys. Lett. 29, 507 (1995); 
J. Magn. Magn. 154, L155 (1996). 

[47] E. P. Raposo and M. D. Coutinho-Filho, Phys. Rev. B 57, 3495 (1998). 

[48] P. Barbosa, E. P. Raposo and M. D. Coutinho-Filho, J. Appl. Phys. 87, 6531 (2000). 

[49] Y. Imry and S. K. Ma, Phys. Rev. Lett. 35, 1399 (1975). 

[50] K. Binder, Z. Phys. B 50, 343 (1983); Phys. Rev. B 29, 5184 (1984). 

[51] D. P. Belanger, A. R. King and V. Jaccarino, Phys. Rev. Lett. 54, 577 (1985). 

[52] A. R. King, V. Jaccarino, M. Motokawa, K. Sugiyama and M. Date, J. Appl. Phys. 57, 3297 
(1985). 

[53] W. Kleemann , Ch. Jakobs , Ch. Binek and D. P. Belanger, J. Magn. Magn. Mater. 177, 209 
(1998). 

[54] D. P. Belanger, A. R. King and V. Jaccarino, Phys. Rev. B 31, 4538 (1985). 

[55] M. Staats, U. Nowak and K. D. Usadel, J. Magn. Magn. Mat. 177-181, 85 (1998). 

[56] J. L. Cardy, Phy. Rev. B 29, 505 (1984). 

[57] S. Fishman and A. Aharony, J. Phys. C 12, L729 (1979). 

[58] H. Rohrer, J. Appl. Phys. 52, 1708 (1981). 

[59] D. P. Belanger, A. R. King and V. Jaccarino, J. Appl. Phys. 53, 2702 (1982); Phys. Rev. Lett. 
48, 1050 (1982). 

[60] D. P. Belanger, A. R. King, V. Jaccarino and J. L. Cardy, Phys. Rev. B 28, 2522 (1983). 

[61] H. Yoshizawa, R. A. Cowley, G. Shirane, R. J. Birgeneau, H. J. Guggenheim and H. Ikeda, 
Phys. Rev. Lett. 48, 438 (1982). 

[62] U. A. Leitao and W. Kleemann, Phys. Rev. B 29, 505 (1984); Phys. Rev. B 35, 8696 (1987). 

[63] A. Aharony, Europhys. Lett. 1, 617 (1986). 

[64] I. B. Ferreira, A. R. King, V. Jaccarino, Phys. Rev. B 43, 10797 (1991); J. Appl. Phys. 69, 
5246 (1991). 

[65] R. J. Birgeneau, Q. Feng, Q. J. Harris, J. P. Hill, A. P. Ramirez, T. R. Thurston, Phys. Rev. 
Lett. 75, 1198 (1995); Phys. Rev. Lett. 77, 2342 (1996). 

[66] D. P. Belanger, W. Kleemann and F. C. Montenegro , Phys. Rev. Lett. 77, 2341 (1996). 

[67] A. R. King, V. Jaccarino, D. P. Belanger and S. M. Rezende, Phys. Rev. B 32, 503 (1985). 

[68] D. P. Belanger, J. Wang, Z. Slanic, S-J. Han, R. M. Nicklow, M. Lui, C. A. Ramos and D. 
Lederman, Phys. Rev. B 54, 3420 (1996). 

[69] M. Lederman, J. V. Selinger, R. Bruinsnia, J. Haniman and R. Orbach, Phys. Rev. Lett. 68, 
2086 (1992). 

[70] U. A. Leitao, W. Kleemann and I. B. Ferreira, Phys. Rev. B 38, 4765 (1988). 

11 



[71] P. PoUak, W. Klccmann and D. P. Belanger, Phys. Rev. B 38, 4773 (1988). 

[72] J. Mattsson, C. Djurburg and P. Nordblad, Phys. Rev. B 61, 11274 (2000). 

[73] U. Nowak and K. D. Usadel, Phys. Rev. B 44, 7426 (1991); Phys. Rev. B 46, 8329 (1992); 
Phys. Rev. B 43, 851 (1991). 

[74] S-J. Han and D. P. Belanger, Phys. Rev. B 46, 2926 (1992). 

[75] S-J. Han, D. P. Belanger, W. Kleemann and U. Nowak, Phys. Rev. B 45, 9728 (1992). 

[76] A. R. King, V. Jaccarino, T. Sakakibara, M. Motokawa and M. Date, J. Magn. Magn. Mat. 
31-34, 1119 (1983); Phys. Rev. Lett. 47, 117 (1981). 

[77] F. C. Montenegro, A. R. King, V. Jaccarino, S.-J. Han and D. P. Belanger, Phys. Rev. B 44, 
2155 (1991). 

[78] C. M. Soukoulis, G. S. Grest, C. Ro and K. Levin, J. Appl. Phys. 57, 3300 (1985). 

[79] J. R. L. de Almeida and R. Bruinsma, Phys. Rev. B 35, 7267 (1987). 

[80] D. P. Belanger, Wm. Murray, F. C. Montenegro, A. R. King and V. Jaccarino, Phys. Rev. B 
44, 2161 (1991). 

[81] K. Jonason, P. Nordblad and F.C. Montenegro, unpublished. 

[82] F. C. Montenegro, K. A. Lima, M. S. Torikachvih and A. H. Lacerda, J. Magn. Magn. Mater. 
177 - 181, 145 (1998); Ising Centennial Colloquium. 

[83] F. C. Montenegro, K. A. Lima, M. S. Torikachvih and A. H. Lacerda, in: Proceedings of the 
Fourth Latin American Workshop on Magnetism, Magnetic Materials and their Applications, 
P.P. MisseU (Ed.), Materials Science Forum, 302-303, 371 (1999). 

[84] A. Rosales- Rivera, J. M. Ferreira and F. C. Montenegro, Europhys. Lett. 50, 264 (2000); 
ICM2000 proceedings. 

[85] J. Satooka, H. Aruga Katori, A. Tobo and K. Katsumata, Phys. Rev. Lett. 81, 709 (1988). 

[86] A. R. King, J. A. Mydosh and V. Jaccarino, Phys. Rev. Lett. 56, 2525 (1986). 

[87] A. E. Nash, A. R. King and V. Jaccarino, Phys. Rev. B 43, 1272 (1991). 

[88] Ch. Binek, S. Kuttler and W. Kleemann, Phys. Rev. Lett. 75, 2412 (1995). 

[89] M. Gofman, J. Adler, A. Aharony, A. B. Harris and M. Schwartz, Phys. Rev. B 53, 6362 
(1996). 

[90] W. C. Barber and D. P. Belanger, J. Appl. Phys. 87, 7049 (2000). 

[91] W. C. Barber, F. Ye, D. P. Belanger and J. A. Fernandez-Baca, unpublished. 



12 




10 
q/c ' 



15 



20 



Figure 1: A comparison of the logarithm of the scaling functions /{q/n) versus q/n for different 
models (See Eq. 1). The pure cases are from approximate expressions from numerical studies M nOJ. 
The REIM and RFIM are determined from the experiments. Note that the corrections to the MF 
equation are largest below the transition and are very significant for the pure d = 2 and random-field 
d = 3 cases. 



13 



<^ 



T3 






5 



4 



3 - 



2 







— 












M u .u 


- 







! 


- h 




- 


a 


1 ~^ 4.5 

-a 

\ 4 


i A 




- 


■ ■ ■ 


J 


, ^ 3.5 


^^ \ 




- 


^F V ^^P B iM^^' 


3 


1 1 1 1 i 1 








70.4 70.6 70.8 


^ 




T (K) 


^ 


^j^tut^k^^'^'^r^ 1 


=^ 


I i ^ 


- 


■1 4 F P" 7 m IF 
1 Tf » r t^o.gs 0.07 2 


_ 


^ T >. \ \ 


_ 


■ ^ ^ ^^■"^^^iS-.J^^ 


~ 


♦ 7 T ^^^^^^^^^^-^^^^^i^^t-^i^ 


- 


^ 7 T (FC) 



70 



72 



74 



T (K) 



Figure 2: Cm vs. T for Fco.93Zno.07F2 at i7 = 7 T as determined using the birefringence technique. 
The inset shows the FC data. There appears to be a tiny hysteresis very close to the transition, 
perhaps a consequence of random-field activated dynamics. 
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Figure 3: Cm vs. T from Monte Carlo simulations modeled after the Feo.8Zno.2F2 system. The 
simularity with the data is striking, though not all Monte Carlo simulations yield a sharp peak in 
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Figure 4: The H — T phase diagram for FexZni-xF2 measured in pulsed magnetic fields. The 
concentrations for the alphabetic labels are the same as in Fig. 5. 
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Figure 5: Low temperature spin flips and phase boundary for FexZni-xF2 as a function of x. Note 
that the phase boundary behavior is quite different for x < Xe and x > Xe- 
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Figure 6: The H ^ T phase diagram for x = 0.56 showing the upper equiHbrium boundary, the 
phase boundary and a lower equiHbrium boundary. 
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Figure 7: Monte Carlo simulation data for the staggered magnetization versus T for magnetic 
concentrations 0.5, 0.6, 0.7 and 0.8. The ZFC and PC procedures exhibit hysteresis for the lower 
concentrations, which only gets worse for slower runs or larger lattices. No hysteresis is observed for 
the higher concentrations. 
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